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ABSTRACT 

We address the deviations of the scahng relations of elhptical galaxies from the ex- 
pectations based on the virial theorem and homology, including the "tilt" of the "fun- 
damental plane" and the steep decline of density with mass. We show that such tilts 
result from dissipative major mergers once the gas fraction available for dissipation 
declines with progenitor mass, and derive the scaling properties of the progenitors. 
We use hydrodynamical simulations to quantify the effects of major mergers with dif- 
ferent gas fractions on the structural properties of galaxies. The tilts are driven by 
the differential shrinkage of the effective stellar radius as a function of dissipation in 
the merger, while the correlated smaller enhancements in internal velocity and stel- 
lar mass keep the slope of the velocity-stellar mass relation near V oc The 
progenitors match a straightforward model of disc formation in ACDM haloes. Their 
total-to-stellar mass ratio within the effective radius varies as Af/A^ cx A^~(0-2-0-3)^ 
consistent with the effect of supernova feedback. They are nearly homologous in the 
sense that the dark-matter fractions within the effective and virial radii scale with 
mass in a similar way, with only a weak decline of density with mass. The progenitor 
radius is roughly R cx A^" '^, compatible with today's intermediate late-type galax- 
ies. This may indicate that the latest dissipative mergers in the history of ellipticals 
involved relatively big discs. The dissipative gas-to-stellar mass ratio is predicted to 
decline oc similar to the observed trend in today's blue-sequence galaxies. 

Such a trend should be observed in relatively massive gaseous galaxies at z 1 — 4. 
The corresponding "baryonic" relations are consistent with homology and spherical 

1/3 

virial equilibrium, V cc R (x M^!^^ . 

Key words: dark matter — galaxies: ellipticals — galaxies: evolution — galaxies: 
formation — galaxies: haloes — galaxies: mergers 



1 INTRODUCTION 

The global properties of elliptical galaxies show non-trivial 
correlations whose origin is one of the most interesting issues 
in the study of galaxy formation. These correlations may 
tell us whether the ellipticals could have indeed formed by 
mergers, and if so what were the properties of their initial 
progenitors. We address these issues in this paper. 

We appeal here to three observable global properties 
of spheroidal stellar systems: the total luminosity L, the 
"effective" radius R encompassing half the total luminosity 
in projection, and the effective velocity dispersion a.^ The 

^ In the literature, the radius is sometimes replaced by surface 
brightness. The use of "effective" surface brightness / = L/B? 
is equivalent to using R, but not exactly so when / and a refer 
instead to the "central" quantities. Nevertheless, as long as the 



galaxies lie on a "Fundamental Plane" (FP) in the three- 
dimensional parameter spac e defined by the logarithms of 
these quantities (Dressier et alj Il987t IPiorgovski fc Davis! 
llQSTtlBernardi fc et al.Jl2003h . with very little scatter about 
it. The projections of this plane onto the planes defined 
by the logarithms of the axes a-L and R-L are well 
fit by straight lines, n amel y "scaling relation s" , following 
iFaber fc JacksonI (1197611 and iKormendvl (Il977l) . 

As specified in the following section, the fundamen- 
tal plane is "tilted" compared to the plane implied by the 
virial theorem. While part of this tilt can b e attributed 
to sy stematic variations in stellar populations dFaber et alJ 
Il987l) . the other part indicates a breakdown of the homol- 



density profile is close to universal, the central quantities should 
roughly scale with the cfTcctivc quantities. 
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ogy, or self-similarity, of the elliptical family, which can be 
expressed as a systematic increase of dark-matter fraction 
within the effective radius as a function of galaxy mass 
jPadmanabhan et al .120041: lOnorbe et alJ2005ll . Another in- 
dication of a breakdown in homology is provided by the 
scaling relations, which imply a steep decline of mean effec- 
tive density with mass, steeper than the moderate decline 
expected based on the propertie s of dark-matter ha loes in 
the standard ACDM cosmology jBullock et al.ll200ll) . Sup- 
porting evidence for non-homology comes from the fact that 
the surface-brightness profiles of larger el lipticals appear to 
be steeper, with a larger Sersic index iCaon et al.l Il993l : 
iGraham fc Guzmanlliool . We seek understanding of these 
deviations from homology. 

In this paper, we address the structural changes induced 
by dissipative mergers. Our simple analytic analysis is based 
on power-law recipes for how the changes in the structure 
parameters due to mergers depend on the effective gas frac- 
tion in the progenitors. We deduce these recipes from a suite 
of disc merger simulations with a varying degree of gas frac- 
tion. This analysis allows us to connect the observed FP and 
scaling relations of ellipticals with the properties of their 
progenitors via simple equations. We use them to put con- 
straints on the properties of these progenitors, their scaling 
relations and especially the systematic trend of gas fraction 
with mass in them. 

In 321 w6 spell out the observed slopes of the scaling 
relations and fundamental plane, and the implied deviations 
from homology. In S|3|we discuss the possible origin of these 
structural deviations, and motivate the study of the role 
of dissipation in mergers. In we formulate the merger 
process in terms of simple power-law relations and derive 
equations that connect the observed relations for elliptical 
galaxies with the properties of their original progenitors, 
including the variation of gas fraction with mass. In ij^we 
use merger simulations to calibrate the necessary power-law 
approximations for the dependence of the structural changes 
in mergers on the gas fraction. In ^ we use the recipes 
derived from the simulations to solve for the properties of 
the progenitors given the observed structural relations. In 
37|we compare the required progenitor properties to what 
might be expected based on other considerations. In fj^we 
summarize our conclusions. 



2 OBSERVED SCALING RELATIONS 

We write the observed scaling relations for elliptical galaxies 
as 



RocL" 



(1) 



where L is the absolute luminosi ty in a given band. The 
slopes d and (3 were determined bv lBernardi et al] i2003l) for 
a large sample from the Sloan Digital Sky Survey (SDSS), 
applying a maximum-likelihood analysis for the conditional 
average of the log of the relevant quantity at a given lumi- 
nosity, namely 



(log crjL) — a log L + const., etc. 



(2) 



Unless specified otherwise, we assume hereafter analogous 
definitions of slopes for all the power-law relations. The most 
recent values based on revised photometry in the r band (M. 



Bernard!, private communication) are a = 0.230 ±0.012 and 
(3 — 0.704 ± 0.025,^ with conditional standard deviations of 
the logs 0.063 and 0.115 respectively. The slopes in the other 
SDSS bands are not very different. 

The Fundamental Plan e as determined from SDSS 
by Bernardi & ot al., (2003) is similar to the original FP 
tDressler et aL 1_9_87. ; .Diorgovski &: Davijll98'^ and can be 
approximated by 



2 r> rH-T 

a Ron L 



(3) 



with f ~ 0.17 and a very small scatter about it. Given the 
way the slopes are determined, the value of f is expected to 
be related to the slopes of the scaling relations via 1 -f f = 
2a -f p. The revised slopes of the scaling relations indeed 
imply T ~ 0.164. 

For the purpose of studying structural and dynamical 
issues we rather replace the luminosity L with the corre- 
sponding stellar mass M (which for ellipticals roughly con- 
stitutes the whole baryonic mass) . We parameterize the sys- 
tematic variation in stellar mass-to-light ratio due to stellar 
populations by 

:^ cx , (4) 

and for the sake of simplicity ignore the scatter about this 
relation. The structural relations for ellipticals within the 
effective radius become 



with 



1 + r- 



(5) 



(6) 



The value of t« is believed to be on the order of 0.1 (e.g. 
Shen et al1 l2003f) . It can therefore make only a ~ 10% dif- 
ference to the scaling slopes a and (3, but a major difference 
to the FP tilt r. 

The scaling relations are non-trivial in two ways. First, 
the FP is tilted relative to the virial plane, 



aR oc M 



(7) 



which is expected to be roughly valid for the total mass 
within the effective radius if the ellipticals are a homologous 
family. The tilt reflects a systematic variation of the total 
mass-to-light ratio within the effective radius, M/L cx . 
This could be partly due to stellar populations and partly a 
structural tilt, 



M 



r = 2q + /3 - 1 . 



(8) 



If indeed r > 0, it represents a breakdown of homology 
for these systems, indicating that the mass fraction of dark 
matter within the effective radius is increasing with mass. 
Furthermore, any non-negative tilt (r ^ 0) for the merger 
remnants is non-trivial because the progenitors are actually 
expected to show a negative tilt relative to the virial plane. 
This is if the variation of M/M with progenitor mass is 
attributed to the decreasing effectiveness of supernova feed- 
back in su ppressin g star form ation as a function of mass 
jDekel fc Silki.l98a-lDekel fc Woa.2003il . 

^ revised from the earlier published values of 0.255 and 0.632. 
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Second, the virial relation also implies that the mean 
total mass density within the effective radius varies with A4 
as 

p oc <x M\ S = 2a ~ 2/3 . (9) 

The observed values of a and (3 correspond to (5 ~ 
—0.95/(1 + r,), namely a rather steep decline of density 
with mass, which provides an even stronger evidence for a 
homology breakdown. The simplest self-similar model of si- 
multaneous spherical collapse predicts virial equilibrium at 
a constant mean density, with a — P = 1/3. Indeed, numer- 
ical studies of dark-matter haloes in ACDM cosmological 
simulations reveal only a weak systematic densit y vari ation, 
weaker than 6 ~ -0.25 (based on Bullock ct al] l200j) . The 
steep decline observed for ellipticals thus implies a break- 
down of homology. Unlike the tilt of the FP, the density de- 
cline cannot be explained by a variation in the stellar M /L 
for any sensible value of r*. 

In our analysis below we adopt r* = 0.1, but this par- 
ticular choice within the range of uncertainty does not have 
a significant effect on our conclusions. The targets for our 
modeling thus become a ~ 0.21 and /3 ~ 0.64, with the 
associated r ~ 0.06 and 5 =~ —0.86. 



3 POSSIBLE ORIGIN OF THE TILTS 

When addressing the possible physical origin of the struc- 
tural tilts, we adopt the standard paradigm of galaxy for- 
mation. It assumes that the accretion of cold gas in dark- 
matter haloes leads to central gaseous discy systems in 
which stars form ( W hitc^fc Recs 1978; Fall & Efstat hioui 
Il980l : IMo et aD Il998h . Major mergers between such sys- 
tems of comparable masses produce more stars, consume 
the available cold gas, and transform the stellar systems 
into elliptical configurations (Fall 1979; J ones fc Efstathioui 



|1979»: Aarseth fc Fall ,1980; Barnes 1992). In haloes more 
massive than ~ IO^'^Mq, there is a shutdown of gas cooling 



^^^^ ^^^^^^ — ^^^^..^^^Q ^^^^ ^ — ^^^^ 

passively and turn "red and dead" (e.g. 


Birnboim & Deke 


12003'; 'Binnov 2004; Kcrcs ct af. 2005; 


Dokcl & Birnboim 


boos-, .Croton et al. 2005; Cattaneo ct al. 


MM- These el- 



lipticals may grow further by non-dissipative mergers with 
other ellipticals. One of the open questions is whether the 
dissipative mergers in the history of big ellipticals involved 
big systems in a relatively late epoch or small systems earlier 
on, followed by major non-dissipative mergers. 

As a simple and general model for the mergers that 
transform the initial systems into ellipticals, we assume that 
progenitor galaxies of type "d" are transformed by mergers 
to remnant galaxies of type "e" . Each galaxy is characterized 
by the three global quantities: a stellar mass M, an effec- 
tive stellar radius R, and a characteristic velocity V (circu- 
lar velocity for discs or velocity dispersion for spheroids, up 
to a multiplicative factor of order unity). The dimension- 
less growth factors describing the structural changes due to 
mergers, Rc/Rd, Vc/Vd and A4e/Md, are determined by the 
merger dynamics. If the progenitors were completely homol- 
ogous, and all the processes involved in a merger were scale 
free, then these factors were independent of mass, and the 
mergers would have led from power-law scaling relations to 



similar power laws, though in general with different zero- 
points. This would have meant that the structural plane in 
M-R-V and its projections are preserved under mergers, but 
perhaps with an increased scatter. 

In particular, if the elliptical galaxies are homolo- 
gous, and the merger dynamics involves gravitational dy- 
namics only, then the FP should be preserved under a 
few major mergers."^ Such a behavior has indeed been 
noticed i n N-body sim ulations of di ssipation l ess m ergers 
iBovlan-K olchin et al.l ^005; Robcrts onet all Hooi). The 
question is thus how did the elliptical galaxies make it to 
the FP m the first place. This requires a breakdown of the 
self-similarity of at least one property of the progenitors or 
the processes acting during the mergers. The dimensionless 
growth factors should depend on mass due to some impor- 
tant merger characteristic that varies strongly with mass. 

The gravitational process associated with spherical col- 
lapse and virial equilibrium are not expected to break the 
self-similarity of the progenitors in a major way. The sys- 
tematic variation of the characteristic collapse time with 
mass introduces a relatively weak variation in halo con- 
centration which is not enough for explaining the devia- 
tions from homology. The origin of angul ar momentum via 
tidal torq ues is also large ly self-similar llFall fc Efstathioui 
[1980; Bullock et alJ 1200 ill , implying disc sizes which are 
close to self-similar. The distribution of merger orbit param- 
eters is not expected to vary significantly with mass either 
JVitvitska et aLll2002ll . 

On the other hand, gas dissipation may lead to a more 
substantial breakdown of homology. Once the progenitors 
contain gas, the dissipation involved in the merger leads 
to a more centrally concentrated stellar system than in the 
non-dissipative case. The tidal interactions during the close 
passages induce collisions and shocks which condense and 
heat the gas. This is followed by enhanced radiative cooling 
and star formation, ending up with a more compact stellar 
system deeper in the potential well. 

Dissipation is indicated, for example, b y the high phase- 
space density in the center s of ellipticals (iHernauist et alJ 
ll993l:lRobertson et al ]200i), and by the tendency for oblate- 
ness and kinematic misalign ment of elliptica ls compared to 
simulated merger remnants JCox et alJl2003) . 

A property of the gas which depends on galaxy mass 
and could therefore produce structural tilts in the scaling 
relations is the cooling efficiency on a dynamical timescale, 
which declines with halo mass (e.g., iRees fc Ostrikedll977l: 
iBlumenthal et al ]|1984V A tilt of this sort has been seen 
in remnants of hydro dynamical merger simulations by 
iRobertson et al.1 ll2006h . Their Fig. 11 shows a tilt in the 
R-M relation when the gas fraction is 0.4 (or higher), but 
the change from /3 ~ 0.44 to 0.51 is small compared to the 
required tilt. A careful inspection of this figure reveals that 
for a fixed gas fraction the tilt is noticeable only for giant 
ellipticals with stellar masses above ~ 5 x lO^^i\f0, while 
the observed tilts extend far below this mass. 



^ Phase mixing and other gravitational relaxation processes may 
eventually push such systems toward the virial plane obeyed by 
the total mass, but this is apparently a long term process that re- 
quires many successive mergers, and, anyway, it does not produce 
a tilt toward the FP. 
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A more pronounced violation of the self-similarity of 
the dissipative processes may arise from a systematic varia- 
tion of the gas fraction with progenitor mass. The relative 
amount of gas in the merging systems has a strong effect 
on the resulting stellar mass and structure of the remnant. 
More gas naturally results in more new stars, and the higher 
degree of dissipation involved in the merger leads to an en- 
hanced shrinkage in the radius R compared to the non- 
dissipative case. Based on a combination of physical con- 
siderations and numerical simulations we have learned that 
the fraction of new stars and the radiative energy loss dur- 
ing the merg er are both roughly pro portional to the effective 
gas fraction iCovington et al]l2006h . These two effects make 
the decline in R/M due to mergers depend strongly on the 
gas fraction. The associated increase in velocity dispersion 
makes the ratio V/M rather insensitive to the mergers and 
the gas fraction. Together, the increase in V is not enough 
to compensate for the strong decrease in R/M, yielding a 
strong gas-dependent decline in the combination relevant to 
the FP, V^R/M. Thus, a proper decline of gas fraction with 
progenitor mass could, in principle, lead to the required tilts. 

Such a gradient of gas fraction with progenitor mass 
is likely for several reasons. We first note that today's disc 
galaxies along the blue sequence show a strong trend which 
can be crudely fitted by a power law, g oc M~'' , where 
q is the gas-to-star mass ratio and 7 ~ 0.4 iKannappanI 
I2OO4I . Figure 2c). The earlier progenitors of elliptical galax- 
ies are expected in general to contain more gas, before it 
has been consumed into stars or removed from the systems. 
Since small-mass galaxies tend to merge earlier than mas- 
sive galaxies in any given neighborhood, the major mergers 
involving small progenitors are expected to be more dissi- 
pative than the mergers of massive progenitors. The gas- 
fraction gradient may also be a result of the higher efficiency 
of quiescent star formation in l arger discs, bein g determined 
by the higher surface density teennicutdllOSSD rather than 
the three-dimensional density. If the mean 3D density is 
roughly independent of mass at any given epoch, as implied 
by the naive spherical collapse model and tidal torque the- 
ory, the surface density is expected to be proportional to 
radius. Finally, a lower star formation rate at lower masses 
is expected to result from the more efficient supernova feed- 
back slowing down the q uiescent star formation preferen- 
tially in smaUer galaxies JPekel fc Silklll98^: iDekel fc Wool 
l2003l:lDekel fc Birnboimll2005^ . 

Hereafter, we take g to refer to the gas that is effectively 
available for dissipation during the merger, i.e., combining 
the initial amount of gas and the dependence of the cooling 
efficiency on mass. 



4 STRUCTURAL TILTS BY MERGERS 

We assume that the progenitors ( "d" ) and merger remnants 
( "e" ) lie on (tilted) planes and scaling relations as in eq. Q 
in the corresponding parameter spaces V^-Rd-Md and Vi- 
Rc-Mc respectively, where all quantities refer to the spheres 
defined by the effective radii encompassing half the stellar 
mass. The six relations can be expressed in terms of the 
growth factors Rc/Rd, Vc/Vd and A4c/Md'- 



QcMf-'^'' . (11) 



and 

(K/l/d) {R./Rd) 
(Mo/Md)''-^ " ' (Xc/Md)A 

We then parameterize the structural growth factors in 
the mergers as functions of g by the power laws 

^'"^ M-d^"' ■ (^2) 

Finally, we assume that the mean dependence of g on pro- 
genitor mass can be parameterized by, 

g^xMr- (13) 

When substituting these power laws in eq. (1101 and 
eqs. (im . we obtain useful relations between the power in- 
dices involved. 



2v + r ^ {1 + Tc)m - 7 ^(tc - Td) 
"fi-v + mac) = (Qc — ctd) , 
^(_r-|-m/3e) = (/3c-/3d) . 



(14) 
(15) 
(16) 



These three equations are degenerate; only two are indepen- 
dent, with T = 2a + P — 1 for "d" and for "e". Once we know 
the growth- factor indices r, v and m from theory, these equa- 
tions connect the final structural tilts to the properties of 
the progenitors, which are described by three independent 
variables: the gas-fraction gradient index 7, and the slopes 
Ud and Pd (one of which can be replaced by Td or Sd)- For 
given values of the observed parameters Oc and /3c (or Tc), 
the solutions for the progenitors are a one-parameter fam- 
ily. Namely, for any assumed value of 7 there is a unique 
solution for Od, Pd, etc. 



MERGER SIMULATIONS: 
DEPENDENCE 



GAS 



\VdJ \RdJ \Md 



(10) 



In order to determine the dependence of the structural 
growth factors on the gas fraction [eq. 11211 1. we ran two sets 
of simulations of typical mergers of disc galaxies, in which 
the initial conditions were identical except for the gas frac- 
tion in the disc; the initial gas-to-stars mass ratio g ranges 
from ~ 0.1 to 3. The two sets differ only in their bulge-to- 
disc ratio. 

The fi ducia l case is G3-G3 from the suite of merger sim- 
ulations in lCoxl i2004ll . Two identical spiral galaxies are put 
on a parabolic orbit, each consisting of stellar and gaseous 
discs and a stellar bulge, all embedded in a ACDM halo. 
The gravitational and hydrodynamical evolution is followed 
using the entropy-conserving, gravitating. S moothed Parti- 
cle Hydrodynamics (SPH) code GADGET dSoringel et alJ 
[2001; Suringcl & Hcrnquist 2002). Gas cool ing, star forma- 
tion and supernova feedback dSpri ngel 200^ are treated us- 
ing simplified recipes that were calibrated to match observed 
star-formation rates. 

The structure of the fiducial progenitor disc galaxies 
mimics Sab-Sb spirals today, similar to the Milky Way. The 
dark-matter halo virial mass is 116x lO^^M©, represented by 
120,000 particles (of 0.97 x IO'^Mq each). The virial radius 
is 272 kpc, corresponding to a vir ial velocity of 135km s~^. 
The halo density profile is NFW llNavarro et alJ 'l997') with 
a concentration parameter C = 6, and a spin parameter A = 
0.05. In the fiducial case, the baryonic fraction of the total 
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mass is 0.054. The total baryons are divided into a stellar 
disc fraction of 0.66, a gaseous disc fraction of 0.20, and a 
stellar bulge fraction of 0.14. The stellar disc mass of 4.11 x 
1O^°M0 is represented by 50,000 particles (of 0.822 x lO^M© 
each). The surface density profile is exponential, with a scale 
radius iid = 2.85 kpc, and the stability is ensured by a 
Toomre Q parameter of 2. The fiducial gas disc of 1.22 x 
1O^°M0 is represented by 50,000 particles (of 2.44 x 10^ M© 
each) . The gas is distributed in an exponential disc of a scale 
radius 37?d extending out to 5 gas scale radii. The fiducial 
bulge is a sphere containing 0.89 x lO'^" Af© of stars in 20,000 
particles. Its density profile is exponential with scale radius 
0.62 kpc. The smoothing length is h = 100 pc for the gas 
and star particles and 400 pc for the dark-matter particles, 
with the force becoming Newtonian at ^ 2.3/i. 

The two galaxies, which were set on a parabolic orbit, 
merge because of dynamical friction due to their massive 
haloes. One galaxy starts such that its disc lies in the col- 
lision orbital plane, and the other galaxy's disc is inclined 
by 30°, both with a prograde orientation of spin. The ini- 
tial separation between the galaxy centres is 250 kpc, and 
the expected pericentre (assuming point-like galaxies) is at 
13.6 kpc, namely about 5% of the virial radius. The merger 
results in two successive starbursts, one after the first close 
approach, and the other after the second, final coalescence. 
The starbursts occur between 2 and 3 Gyr after the begin- 
ning of the simulation, and the remnant is "observed" at 
4 Gyrs. The amount of stars formed during the merger is 
roughly proportional to the initial gas fraction, and is not 
too sensitive to the orbit or orientation. In the cases with 
moderate gas fractions, the young stars formed during the 
merger typically constitute ^ 20% of the total stars. 

The remnant galaxies resemble nor mal elliptica l galax - 
ies in many respects, as demonstrated in lDekel et al.l i2005ll . 
This includes the stellar density profile, the ellipticity and 
triaxiality, the velocity-dispersion profile as well as higher 
moments of the line-of-sight velocity distribution. 

The galaxies in all cases are identical in their dark mat- 
ter halos, baryonic fraction and merger parameters. In one 
suite of simulations the bulge-to-disc ratio is the fiducial 
0.17, while in the other suite it is 0.39. Each of the two suites 
consists of 5 cases, where the bulge-to-disc ratio is kept fixed 
and only the disc gas-to-stars mass ratio g is varied. In the 
small-bulge suite, g varies between 0.12 and 3.0 [correspond- 
ing to gas/(gas-|-stars) between 0.11 and 0.75], while in the 
suite with a more massive bulge, g varies between 0.097 and 
1.72 [corresponding to gas/(gas-|-stars) between 0.088 and 
0.63]. 

Figure Q shows the growth factors in the mergers as a 
function of g. The trends are fit by the power laws in eq. 11211 
with r ~ —0.4, w ~ 0.12, m ~ 0.3. The radius growth factor 
shows the strongest dependence on g. The power laws are 
good approximations for the g-dependence of the radius and 
velocity factors over the whole range of g, but the power-law 
fit in the case of stellar mass is rather crude. The effective 
slope of m ~ 0.3 coincides with the local slope near g ~ 0.5, 
and it can serve as a sensible approximation over a g-range of 
somewhat less than a decade about it. At very low g values 
the slope is as small as ~ 0.1, and at very high g values it 
becomes steeper than 0.5. Still, the crude power-law fit can 
adequately serve our simplified modeling of the tilts. 

While the magnitude of the radius shrinkage depends 




0.4 : 




-0.2 V . 

-1.0 -0.5 0.0 0.5 

Log g (gas/stars) 

Figure 1. Dependence on gas content in merger simulations. Tiie 
relative changes in effective radius R, internal velocity V and 
stellar mass M are plotted against the initial gas-to-star ratio 
g The two curves, which refer to different bulge fractions, show 
similar slopes. The effective power-law fits are indicated. 

on the bulge-to-disc ratio, we find that the power-law fits 
are insensitive to the presence of a more massive bulge. This 
indicates that our results should be robust to the exact shape 
of the progenitors. 

When the growth factors are plotted instead as func- 
tions of g' =gas/(gas-|-stars) rather than g, the correspond- 
ing best-fit power indices become r' ~ —0.6, v' ~ 0.2, 
m' ~ 0.4. However, the power-law fits in this case are not as 
good. In particular, the R factor flattens at the low-i;' end, 
and the m factor steepens significantly at the high-g' end. 
We therefore prefer to use g as the gas-content variable. 



6 RESULTS: PROGENITOR PROPERTIES 

We are now set to solve eqs. 1141 - 116j . We use the values of 
the growth-factor indices r — —0.4, v = 0.12 and m = 0.3 
as deduced from the merger simulations. We then adopt 
the observed values for ellipticals So = 0.23 and /3o = 0.70 
(fc = 0.16) and solve for the required progenitor properties. 
The one-parameter family of solutions is shown in Fig.|5| For 
each choice of 7, the associated values of the power indices 
Qd and Pd, and their linear combinations Td and Sd, can be 
read from the corresponding curves. For example, the verti- 
cal line at 7 = refers to the null case where the progenitors 
are identical to the final ellipticals because the available gas 
fraction does not vary with mass. Any other vertical line 
corresponds to another possible solution for the progenitor 
properties. Higher initial values of 7 naturally correspond to 
larger changes in the power indices as a result of mergers. 
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Figure 2. Progenitor properties (the subscript "d" is omitted) 
that load by mergers to the observed structural indices of ellip- 
ticals, «e = 0.23 and = 0.70 (j^ = 0.16). The values of the 
parameters along each vertical line correspond to a solution of 
eqs. IT41-IT61. The power index 7 refers to the decline of gas- 
fraction with mass. The indices a and fi are the slopes of the 
V M. and KM. scaling relations. The tilt r refers to variations in 
the ratio of mass to stellar mass within the effective radius. The 
index (5 refers to variations in mean total density within the effec- 
tive radius. The solid and dashed curves correspond to r« = and 
T« = fe, where the tilt of the FP compared to the virial plane is 
either all structural or all due to stellar populations respectively. 

We see that the slope a of the V M relation does not 
change much by dissipative mergers even if 7 is large. This 
is because both V and M increase during the merger in a 
way that makes the systems propagate roughly along the 
V M relation. Also, given that v = 0.12 is so small, and that 
both 7 and m are most likely smaller than unity, eq. I|15|l 
says that is just a bit larger than Ofc, quite insensitive 
to the actual value of 7. The observed slope of Qc — 0.23 
thus puts a strong robust constraint on the corresponding 
progenitor slope, forcing it to be well within the range 0.20- 
0.28. This rules out in particular progenitors that lie on the 
naive virial-spherical plane ad = /9d = 1/3 (but we will see 
below that this conclusion is limited to the scaling relations 
based on stellar mass, and becomes invalid when all the 
baryonic mass is considered). 

On the other hand, the change in the slope /3 of the 
RM relation is strongly dependent on 7 and can be quite 
large. This is because the radius shrinkage is the dominant 
change, helped by the fact that r and m in eq. (Ilfcil are of 
opposite signs. Correspondingly, the values of r and 5 suffer 
changes which are strong functions of 7: the former grows 
and the latter declines as a result of mergers. 

If we limit the discussion to conventional models of 
galaxy formation which do not permit an increasing den- 
sity with progenitor mass, the corresponding requirement 
5d ^ implies an upper limit of 7 ^ 0.67 (for t* = 0.1). 
This corresponds to /3d ^ 0.25 and rd ^ —0.25. 

If we also require that the initial M/M is not increasing 



with progenitor mass, based for example on the common 
wisdom that su pernova feedback b ecomes less effective in 
massive haloes (iDekel fc Sil3ll986l) . then the requirement 
T"d ^ puts a lower bound of 7 ^ 0.15 ((for r« = 0.1). This 
corresponds to /3d ^ 0.55 and 5d ^ —0.7. 

Progenitors with 7 ~ 0.4 — 0.6 are favorable not only 
because these values are compatible with today's trend of 
gas fraction along the blue sequence, but also because they 
represent solutions which are close to homology, roughly 
obeying the virial theorem (rd ^ 0) and a weak decline of 
density with mass (5d ^ 0). For example, a typical value 
of 7 = 0.55, compatible with today's gas gradient along 
the blue sequence plus an estimated contribution from the 
mass dependence of the cooling efficiency, corresponds to 
ad ~= 0.25, /?d =^ 0.31, Td ~ -0.20 and 5^ ~ -0.15. 

We learn from Fig. |21 that the value of r, , namely the 
way the stellar M/L varies with mass, has a surprisingly 
small effect on the required progenitors. Even a choice of a 
large stellar tilt, t« = fe = 0.16, which eliminates the need 
for a structural tilt of the FP (namely To = 0) , still implies 
strong changes in /3 (and r and 5) between the progenitors 
and the remnants once 7 is not negligible. 



7 INTERPRETING THE PROGENITORS 

The solutions obtained above can be compared to what 
might be expected based on other considerations for the 
properties of the early gaseous progenitors of ellipticals. 

Recall that the naive model of top-hat collapse to virial 
equilibrium followed by disc formation conserving angular 
momentum, at the same epoch for all galaxies, yields ad = 
/?d = 1/3, namely progenitors that lie on the virial plane 
(rd = 0) with the same effective density for galaxies of all 
masses ((5d = 0). As noticed in Fig. |21 this is not a proper 
solution. A solution with /3d = Od does exist, for 7 ~ 0.67, 
but it has the "wrong" value of /3d = ad — 0.25. A value 
close to ad — 1/3 cannot be obtained for any sensible value 
of 7 as long as the observed slope is ^ 0.23, because a 
hardly changes in the course of a merger. 

However, a somewhat more sophisticated conventional 
model of disc formation in ACDM haloes does predict struc- 
tural power indices in the range required for the progen- 
itors of ellipticals. Studies of the profiles of d ark haloes 
in cosmological si mulation s jBullock etaPbOOj) . based on 
the NFW (Navarr o et al.lll99 7) functional form (but with- 
out limiting the generality of the analysis), yield a mean 
trend of halo concentration with virial mass, Cv oc M~^'^''^ , 
reflecting the diff erent effective formation epoch as a func - 
tion of halo mass llBullock et alj200ll:IWechsler et al.l2002l) . 

While the flat part of the circular velocity resembles the 

1/3 

virial velocity, Vv oc Mv , the maximum NFW circular ve- 
locity is Vmax oc c^'^'^Vv (X A^v '^''. If the relevant velocity 
characterizing the progenitors is between 14 and Vmax, the 
model predicts a slope for the Vd-A^d relation in the range 
0.30 ^ ad/(l -I- Tv) ^ 1/3, where Tv characterizes possible 
variations of virial to stellar mass ratio, My/ M oc A^^". 

T he effective disc radius is evaluated bv lBuUock et all 
i200ll) based on the dark-halo properties under the assump- 
tions of a constant spin parameter as a function of mass (as 
seen in the simulations) and conservation of angular momen- 
tum during the gas contraction. For a disc baryonic fraction 
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that is independent of mass, and in the hmit Cv 3> 1, the 
disc radius can be approximated by J? oc c^^'^-Rv oc M^''^^ . 
An increasing baryonic fraction with mass would tend to 
counter-balance the c7°'^ dependence on concentration, so 
the model predicts a slope for the -Rd-A^d relation in the 
range 1/3 13^/(1 +t^) «? 0.42. 

This crude model based on the ACDM halo properties 
is compatible with solutions shown in Fig.|5|for the required 
properties of the progenitors of ellipticals. If we take the 
model predictions to be ad — 0.32(1 -|-Tv) and /3d — 0.37(1 + 
Tv), and assume as before r* ~ 0.1, then the proper solution 
of eqs. Jinj andiniis Tv ~ -0.24, 7 ~ 0.6, ad = 0.24 and 
/3d = 0.28. The corresponding tilt of M/A4 is negative, Td — 
—0.23, and the density decline is rather weak, Sd — —0.08. 

The value of 7 ~ 0.6 is roughly compatible with today's 
gas gr adient along the blue sequence of 7 ~ 0.4 lKannappanl 
i2004l) . The difference can be attributed to the mass depen- 
dence of the cooling efficiency, when g represents the gas 
that actually manages to cool on a dynamical time scale 
after the merger. This contribution to 7 can be crudely es- 
timated from the steepening from /3d — 0.52 to /3e — 0.44 
seen in Fig. 11 of iRobertson et aP l|200^ for mergers with 
a fixed gas fraction g' = 0.4 (corresponding to 17 = 0.67). 
When inserted in eq. (I16II . with r — —0.4 and m — 0.3, we 
obtain a crude estimate of 7 ~ 0.15 for the cooling-efficiency 
contribution at a fixed g. A value of 7 ~ 0.5 — 0.6 for the 
gaseous progenitors is therefore not surprising. 

Unlike the final ellipticals, the gaseous progenitors are 
predicted to be self-similar in the sense that the ratio of 
total-to-stellar mass scales with mass in a similar way inside 
the virial radius and inside the effective radius, Tv — Td. 
The fact that this tilt is negative, ~ —0.24 is not surprising 
either. Such a tilt is seen below ~ 3 X 10^° Mq iBell et alJ 
l2003l:lYang et all2003f) . A negative tilt as steep as Tv ~ —0.4 
is deduced from the "fundamental line" of dwarf galaxies 
and predicted by crude energy considerations of supernova 
feedback in the limit where a small fra ction of the gas has 
turned into stars very eff ective feedback jPekel fc Woo 2 0o3 
iDekel fc Birnboimll2005l) . 

The slope /3d — 0.28 lies well inside the range of the 
slopes observed for late-type galaxies in the SDSS: /3d — 0.15 
and 0.40 below and above ~ 4 x W^'^Mq respectively 
(Shcn ct al. 2003; Kauffmarm et al. 2003). A similar slope 
of /3d — 0.29 is also deduced from a study of a large local 
sample o f late- type spirals, ado pting M/Lj oc Lj'^^ for these 
galaxies (ICourteau et al.ll2006^ . The R-Ai relation for typi- 
cal progenitors thus resembles that of today's intermediate, 
~ I/, spirals. 

The value of Qd — 0.24 is compati ble with the slop e 
of the "baryonic" TuUy-Fisher relation bv lMcGauehl J2005ll . 
referring to the flat part of the rotation curve. It is perhaps 
a bit lower than the Od — 0. 28 deduced using M/Li oc L?'^^ 



from ICourteau et al.l l|200^, in which the rotation velocity 
is measured at 2.2 exponential disc radii. 
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Figure 3. Same as Fig.^but for the baryons (gas+stars) replac- 
ing the stars. The baryonic mass is preserved during the mergers. 
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Figure 4. Same as Fig.|2]but for the baryons (gas+stars) replac- 
ing the stars, showing the solutions to eqs. 1171 . 



the fact that the change in Mbar during the merger can be 
neglected, namely m' = in the analog of eq. (1121 . It also 
eliminates the need to deal with the poor power-law fit of 
Mc/Md as a function of g (Fig.0. The analogs of eq. 1151 
and eq. 11161 become the explicit expressions 



Qd = fto + " 7 



/3d = /3c + r 7 



(17) 



8 BARYONIC RELATIONS 

The analysis above could be performed alternatively with 
the stellar mass M replaced by the baryonic mass Mbar (of 
stars plus gas). We use the same notation as above with a 
prime added to each variable. This analysis is simplified by 



in which 7' refers to the gas-to-baryon ratio in the progeni- 
tors, 



g' oc Af, 



7 

bar 



(18) 



Figure |H] is the analog of Fig. Q now as a function of 
the gas fraction g' , showing reasonable power-law fits with 
slopes v' = 0.17 and r' = 0.42. 
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The solutions of eqs. 11711 are shown in Fig. |1| With 
m' = 0, the changes in a' and S' as a function of 7' are 
stronger than they were in a and 5 as a function of 7, while 
the changes in /3' and r' are weaker than before. We see in 
particular that a proper solution is ~ ~ 0.33, namely 
— 5d — Oi with 7' ~ 0.7. This is the only solution for 
which neither the effective total-to-baryonic mass nor the 
effective density are increasing with mass. If the baryonic 
mass scales with the virial mass, these progenitors are the 
homologous family predicted by the simplest model of spher- 
ical collapse into virial equilibrium. This indicates no signif- 
icant loss of baryons from the progenitors due to feedback 
or other effects. The feedback may have an important role 
in determining the strong variation of gas-to-baryon ratio 
with progenitor mass. 



9 CONCLUSION 

We have clarified how the deviations of the structural scaling 
relations of elliptical galaxies from the expectations based on 
the virial theorem and an assumed homology of the elliptical 
family could arise from the dissipative processes in major 
mergers. The key is a systematic variation of the gas fraction 
available for dissipation as a function of progenitor mass. 
The tilt of the elliptical's fundamental plane compared to 
the virial plane and the steep decline of density with mass 
are driven by the differential shrinkage of the effective stellar 
radius as a function of the degree of dissipation involved in 
the merger. 

The correlated enhancements in internal velocity and 
stellar mass during a merger, and the fact that the change 
in velocity is rather small, lead to only a small change (~ 
10%) in the log-slope of the velocity-stellar mass relation, 
keeping it near V (x Ai^^'^ for both the merger progenitors 
and products. This is in the ball-park of the ste llar TuUy- 
Fisher relation observed in today's disc galaxies iMcGauehl 
I2OO5I: ICourteau et aljlioo^ . The small changes in velocity 
indicate that the tilts are primarily associated with structure 
rather than kinematics. 

The predicted log-slope of the radius-stellar mass re- 
lation in the progenitors is roughly R oc M"'^, compatible 
with the observed relation for intermediate late-type galax- 
ies ne ar M - 10^°" "Mg JShen et alj|2003l: ICourteau et alJ 
|200(J) . It is steeper than the slopes observed at the faint end 
{f3d — 0.15) and at the bright end (/3d — 0.4) of late-type 
galaxies. 

The similarity of the progenitors to today's relatively 
big discs may imply that the latest dissipative mergers in 
the history of ellipticals tended to involve grown discs rather 
then smaller dwarfs. This is compa tible with the detection s 
of big gaseous discs at 2 ~ 2-4 fTacconi & et al.,"2006"), 
which could be produced by cold flows ( Dekel fc Bi rnboim 
[200^, and with the histories of big ellipticals in cosmolog- 
ical simulations that prope rly match the galaxy bimodality 
at low and high redshifts dCattaneo et alj|2006l . Cattaneo, 
Dekel & Faber, in preparation). However, it does not rule 
out later non-dissipative mergers along the red sequence, 
which tend to preserve the FP and the associated scaling 



relations.'' The tilts can be explained by dissipative major 
mergers, or related dissipative processes, somewhere in the 
histories of elliptical galaxies. 

The required progenitor properties match the expecta- 
tions from a simple model of disc form ation in ACDM haloes 
JMo et aklllioSlBullock et al.ll200J) . First, the progenitors 
are predicted to he on a plane that is tilted relative to the 
virial plane in the opposite sense to the tilt of the fundamen- 
tal plane of ellipticals, roughly M/M oc 7\/(-(0-2-0-3) , This is 
similar to the expectations from the effect of supernova feed- 
back. The tilt becomes positive as a result of the dissipative 
shrinkage of the baryonic component in the mergers. 

Second, The progenitors are expected to be close to ho- 
mologous, in the sense that the dark-matter fractions within 
the effective radius and the virial radius scale similarly with 
mass. They should have a very weak decline of density with 
mass, roughly p cc A^""'', which later steepens by the differ- 
ential dissipative shrinkage to near p oc M-"-^ in the merged 
ellipticals. 

Our main prediction is that the variation with pro- 
genitor mass of the dissipative gas-to-star ratio should be 
roughly g cx M~'' with 7 ~ 0.5-0.6. This is consistent with 
the observed trend (7 ~ 0.4) in today's blue-sequence galax- 
ies iKannaDDanll20041 . after subtracting a certain contribu- 
tion to 7 from the dependence of cooling efficiency on mass 
at a fixed gas fraction. In an ongoing study we attempt to 
disentangle in more detail between the roles of gas-fraction 
and cooling efficiency. This will be achieved by performing 
merger simulations similar to the ones described in this pa- 
per but between galaxies of higher and lower masses. 

The corresponding baryonic relations, in which the stel- 
lar mass is replaced by the total baryonic mass in gas and 

stars, are consistent with the simplest homology and spher- 

1/3 

ical virial equilibrium, V oc R oc M^^i^^ , namely negligi- 
ble variations in M/Mbar and in mean density. While the 
baryonic fraction is hardly changing, the gas fraction is pre- 
dicted to vary strongly with baryonic mass in the progeni- 
tors, g' cx Af-o/. 

The little relevance of the bulge fraction to the slopes of 
the scaling relations indicates that the gaseous progenitors 
need not necessarily be pure discs, as long as they have a 
proper systematic decline of gas fraction with mass. 

Our qualitative conclusions are insensitive to the actual 
contribution of variations in stellar mass-to-light ratio to the 
tilt. A substantial variation in the stellar mass-to- light ratio, 
on the order oi M/L oc Ad'^'^ , could in principle contribute 
to the tilt of the fundamental plane, but it cannot give rise 
to the observed strong decline of density with mass. This 
indicates that the structural changes due to differential dis- 
sipation effects in mergers should have an important role in 
shaping up the properties of ellipticals, and implies that the 
trend in available gas fraction as a function of progenitor 
mass is a robust prediction. 

High-redshift observations should identify the predicted 
early progenitors of today's ellipticals, namely relative mas- 
sive, gas-rich discy galaxies with a significant decline of gas 
fraction as a function of mass. Indeed, an analysis of the mor- 
phology of z ~ 1 galaxies using a non-parametric classifica- 



* This is yet to be tested using simulations of several successive 
non-dissipative mergers between spheroids. 
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tion system jLotz et al. suggests that a large fraction 

of these galaxies are massive discs. Since this class of object 
makes up a significant fraction of the star formation, they 
are likely to be gas-rich. Observations at higher redshifts in 
the range z ~ 2-4, both in the UV and sub-millimeter wave- 
lengths, have revealed gas rich galaxies of ~ 10^^ Lq, whose 
kinematics rese mble in many cases ro tating thick discs of 
radii ~ 2-3kpc iTacconi fc et al!1l2006ll . These galaxies may 
represent the required progenitors. A straightforward test of 
the differential dissipation model would be the detection of 
the predicted systematic decline of gas fraction with mass 
in such high-redshift galaxies. 
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